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algorithm 
0 引言 能 直接 用 以 剩余 寿命 预测 。 总 体 而 言 ， 现 阶段 HMM 模型 的 
尽 研究 因 其 自身 限制 也 仅 停 留 在 设备 健康 诊断 方面 。 
隐 马 尔 可 夫 模型 (hidden Markov model，HMM) 作 为 一 种 为 了 克服 HMM 模型 的 局 限 性 ， 学 者 们 提出 了 隐 半 马尔 
统计 学 方法 ， 从 概率 角度 描述 系统 外 在 表现 与 隐 含 状态 之 间 可 夫 模 型 (hidden semi Markov model HSMM)，Donglsl 等 运用 
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基于 改进 退化 隐 马 尔 可 夫 模 型 的 设备 健康 诊断 与 寿命 预测 研究 
刘 文 溢 ， 刘 勤 明 '， 叶 春明 ， 李 冠 林 


(上 海 理 工大 学 管理 学 院 ， 上 海 200093) 


摘 要 : 针对 隐 马 尔 可 夫 模 型 在 进行 设备 健康 诊断 时 与 实际 存在 较 大 偏差 的 问题 ， 提 出 了 一 种 以 似 短 关系 加 速 退 化 
为 核心 的 改进 退化 隐 马 尔 可 夫 模 型 (DGHMM)。 首 先 ， 引 入 退化 因子 描述 设备 衰退 过 程 , 提出 的 似 回 关系 加 速 退化 较 
常规 指数 式 加 速 退化 而 言 ， 能 更 好 地 描述 设备 服役 期 间 随 着 役 龄 增加 性 能 的 逐步 下 降 。 其 次 ， 以 全 局 搜索 能 力 相 对 
较 强 的 改进 遗传 算法 代替 常规 EM 算法 进行 参数 估计 ， 克 服 了 EM 算法 易 陷 入 局 部 最 优 的 局 限 性 。 同 时 ， 针 对 隐 马 
尔 可 夫 模 型 时 间 上 须 服 从 指数 分 布 而 不 能 直接 用 以 寿命 预测 的 局 限 性 问题 , 提出 了 一 种 以 近似 算法 与 Viterbi 算法 为 
基础 的 贪 禁 近 似 法 ， 以 寻求 最 大 概率 剩余 观测 为 目的 ， 动 态 地 寻求 最 大 概率 剩余 状态 路 径 ， 对 设备 剩余 寿命 进行 预 
测 。 最 后 ， 通 过 美国 卡特 彼 勒 公司 液压 泵 数据 集 对 所 提出 的 方法 进行 验证 评价 。 结 果 表 明 ， 基 于 改进 退化 隐 马 尔 科 
模型 的 设备 健康 诊断 与 寿命 预测 方法 在 描绘 设备 退化 、 设 备 状态 诊断 准确 率 方 面 更 加 有 效 ， 在 剩余 寿命 预测 上 亦 为 
可 行 。 

关键 词 : 隐 马 尔 可 夫 模 型 ,设备 退化 ;健康 诊断 ;剩余 寿命 预测 ; 遗传 算法 ; 近似 算法 
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Equipment health diagnostics and prognostics method based on improved degenerated HMM 


Liu Wenyl Liu Qinming, Ye Chunming, Li Guanlin 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: In order to solve the problem of large deviation between Hidden Markov Model and actual equipment health 
diagnosis, this paper developed an improved Degenerated Hidden Markov Model (DGHMM) with a core of the quasi power 
relation. First, the model adopted the degradation factors, modeling the process of recession for the equipment’s continuous 
decrease in performance. Compared with the conventional exponential accelerated degradation, the quasi power relation 
accelerated degradation can better describe the process that the performance of the equipment decreases gradually with the 
increase of service age. Then, the improved genetic algorithm can replace the conventional EM algorithm for parameters’ 
estimation, which overcomes the limitation that the EM algorithm is easy to fall into local optimization. At the same time, in 
terms of the limitation of life prediction problem as a result of the Hidden Markov Model must obey exponential distribution, 
an algorithm named greed & approximation based on approximation algorithm and Viterbi algorithm came out, and to seek 
maximum probability remaining observation, for the purpose of seeking maximum probability dynamically surplus state path, 
to predict the residual life of equipment. Finally, the proposed method is validated and evaluated with the data set of caterpillar 
hydraulic pumps. The results show that the method of equipment health diagnosis and life prediction based on the improved 
degraded hidden Markov model is more effective in describing equipment’s degeneration and the accuracy of equipment state 
diagnosis, and is also feasible in the prediction of residual life. 

Key words: hidden Markov model; degradation health diagnostics; residual useful lifetime; genetic algorithm; approximation 


的 关系 ， 既 能 反映 对 象 的 随机 性 ， 又 能 反映 其 潜在 结构 ， 适 用 HSMM 的 状态 驻 留 对 设备 进行 了 健康 诊断 与 剩余 寿命 预测 ; 
复杂 系统 的 建 模 问 题 ， 不 仅 在 语音 识别 及 分 词 方面 取得 Liu 运用 K-means 算法 及 交叉 验证 算法 确定 电网 退化 阶段 ， 


YH 


于 处 到 


了 应 用 驻 汪 ， 在 模式 识别 与 健康 诊断 方面 也 取得 了 良好 的 效果 。 ”并 用 HSMM 对 电网 进行 诊断 及 寿命 预测 ， 除 此 之 外 对 


CareyBl 等 通过 分 析 机 械 设 备 的 振动 信号 首次 将 HMM 成 ”多 传感器 ， 多 输入 的 模型 ， 学 者 们 也 提出 了 一 系列 方法 


基于 


整合 


功 应 用 于 故障 诊断 领域 ，Wang 欠 等 建立 了 包含 耦合 的 离散 时 ”传感器 信号 ， 使 得 模型 能 够 更 准确 地 进行 系统 的 状态 识别 及 


变 时 消 和 耦合 的 分 布 时 变 时 滞 的 模型 运用 有 效 技术 约束 不 确 。 寿命 预测 站。 上 述 研究 有 着 一 个 共同 点 ， 即 HSMM 本 
定 项 ， 并 通过 数值 仿真 验证 了 方法 有 效 性 。Mab] 等 阐述 了 赖 于 先 验 分 布 , 从 而 不 可 避免 地 存在 一 定 程度 的 主观 性 特征 。 
HMM 在 时 间 上 存在 必须 强制 符合 指数 分 布 的 不 合理 性 而 不 随 着 对 HMM 及 HSMM 的 研究 逐渐 深入 ， 研 究 者 
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运 而 生 
实际 上 随 3 


a Chen 等 指 


刘 文 溢 ， 


: 基于 改进 退化 隐 马 尔 科 夫 模 型 的 设备 健康 诊断 与 寿命 预测 研 


现 马尔 可 夫 类 模型 对 于 状态 的 识别 与 实际 依然 存在 较 大 的 差 
异 ， 为 了 得 到 更 符合 设备 实际 情况 的 模 
传统 模型 状态 转移 概率 是 固定 的 ,但 


设备 工作 时 间 的 增加 ， 设 备 损耗 逐渐 积累 ， 设 备 


型 ， 设 备 退化 理论 应 


的 状态 转移 必然 发 生变 化 ， PengL 设 计 了 三 类 退化 因子 加 入 


到 HSMM 模型 当中 ， 并 对 加 入 退化 


了 重新 推导 ， 对 参数 进行 了 寻 


用 到 HSMM 模型 
化 的 有 效 性 。 
到 模型 中 并 


上 述 对 


于 设备 退化 的 下 
段 设 不 同 的 先 验 分 布 ，T 
进 优 化 的 研究 还 暂时 空缺 。 

极 大 似 然 估 计 作 为 参数 估计 的 常用 方法 ， 通 过 
验 数 据 出 现 似 然 值 最 大 的 参数 ， 来 对 
参数 进行 合理 估计 。 相 应 地 ，EM 算法 以 其 相对 简单 ， 
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生 的 概率 是 不 变 的 ， 而 这 个 假设 限制 了 HMM 实际 建 模 分 析 
中 的 能 力 。 对 此 , 国内 外 学 者 提出 了 时 变 HMM。 典型 的 设备 
退化 规律 可 由 图 1 所 示 。 整 个 退化 过 程 可 分 为 A、B、C 三 个 
阶段 ， 其 中 A 段 为 平稳 退化 阶段 ， 此 阶段 设备 性 能 指数 基本 
呈现 稳定 态势 ，B 段 为 均匀 退化 阶段 ， 设 备 性 能 指数 呈现 均 


妹子 的 模型 的 算法 进行 


E 估 ; Hen 等 将 三 类 退化 因子 运 
型 , 并 搭建 柱 塞 泵 实验 台 , 通过 实例 验证 了 退 


究 ， 大 多 是 将 退化 应 用 


高 的 特点 而 被 广泛 应 | 


易 陷入 局 部 最 优 的 缺点 。 
将 智能 算法 应 用 到 参数 估计 领域 .Yuan 


于 参数 估计 ， 


Yu 等 将 遗传 算 


语法 运用 至 


对 于 退化 因子 本 身 的 改 


寻找 使 实 
一 些 复杂 模型 中 的 未 知 
效率 
但 EM 算法 本 身 却 具有 


为 了 寻求 拓展 ， 越 来 越 多 的 研究 者 


3 等 人 | 


改进 PSO 算 


法 对 Baum-Welch 算法 进行 修正 ， 获 得 了 更 加 精确 的 参数 ; 
到 灰色 模型 中 进行 参数 估计 , 数值 结 


果 显 示 ， 采 用 
用 基于 


HMM 


GA 可 大 大 提高 模型 预测 精度 
遗传 算法 的 函数 逼近 器 估计 感应 电机 的 参数 ， 并 取得 
了 较 其 他 方法 更 好 的 效果 ; 
模型 进行 优化 ， 得 到 了 辨识 度 最 高 的 初始 观测 矩阵 ; 


Zhanglg 等 运 


; Krause05 等 采 


人 工 免 疫 算法 对 


Zhangl7 等 运用 自 适 应 基因 粒子 群 算法 优化 隐 马 尔 可 夫 模 型 


(CHMM), 得 到 了 比 传统 HSMM 更 精确 的 结果 。 
为 全 局 搜索 能 力 较 强 的 智 


类 模型 中 ， 存 在 着 可 延展 开拓 的 领域 。 


9 能 寻 优 算法 ， 尚 未 应 用 到 马尔 可 夫 


遗传 算法 作 


基于 上 述 问题 


页 ， 首 先 引入 退化 因 


模 ， 得 到 设备 退化 过 程 的 描述 ， 并 提出 了 一 种 以 似 窜 关系 加 


速 退化 为 核心 的 DGHMM， 较 常规 指数 式 加 速 退化 而 言 ， 能 


更 加 准确 地 描述 设备 1 


牛 


以 改进 遗传 算法 代替 常规 卫 
针对 HMM 模型 不 能 
一 种 基于 近似 算法 与 Viterbi 算法 的 贪 禁 


同时 ， 
提出 了 


子 ， 对 设备 衰退 进行 建 


能 随 役 龄 增加 而 逐渐 加 速 下 降 的 过 程 。 
M 算法 对 整个 模型 进行 参数 估计 
直接 用 以 寿命 预测 的 局 限 性 问题 


村 近似 法 》 以 设 


备 退 化 时 状态 转移 概率 矩阵 的 变化 为 载体 去 预测 设备 剩余 寿 


让 


命 。 最 后 ， 


1 HMM 模型 


HMM 模型 的 表达 式 可 描述 为 4 = (x, 4, B) , 模型 参 炎 
a) N : 状态 数 , 状态 所 构成 的 集合 
在 任 一 时 刻 ! 所 处 的 状态 为 9 ， 


以 美国 卡特 彼 勒 公 
证 评价 ， 结 果 显 示 ， 本 文 方法 是 可 行 的 和 


司 液压 泵 数据 集 对 模型 进行 验 


效 的 。 


由 术 如 下 : 
为 {5153,54…Sy 41S} ， 


S, E (Si 92 83" SN 1 SN) 。 


b)M : 观测 数 , 观测 所 构成 的 集合 为 {91,&,&…9wa,0w} ， 


Cc) TX: 


有 概率 之 和 为 1， 初 始 概率 向 量 


任 一 时 刻 ! 状态 i 产生 的 观测 为 0, @ (89,8,@… 
初始 状态 概率 分 布 ， 起 三 (M1, To, TAN Tv)， 其 中 所 
的 意义 为 HMM 开始 时 ， 即 


-Or Or) o 


1 = 工时 ， 模 型 处 于 第 was 山 23…N-LN) 个 状态 的 概率 。 


d) 4 :状态 转移 和 矩阵 。4 是 一 个 Nx 的 矩阵 , A=[& 
状态 7 的 概率 ， 


ww 为 从 状态 转移 到 


中 ww 7 


BN w =PG ,=7j1S = 有 且 


D0 le 23 NIN), D0,=1j ed,2,3N-LN). 


e) B : 观测 


b(K)= Plo, =0,|S,= )(j (1,2,3.N), K e (1,2,3:..M)), 


2 ”设备 退化 过 程 
2.1 退化 因子 设计 


概率 分 布 。 


其 中 
且 > DOD=1 
k=] 


B=[b,(K)] 


NxM 9 


传统 HMM 假定 在 设备 整个 寿命 周期 内 ， 其 状态 转移 发 


匀 变 化 ; C 段 为 加 速 阶段 ， 设 备 性 能 随 着 寿命 的 推移 呈现 加 


速 下 降 的 趋势 ， 健 康 状况 急剧 恶化 ， 直 至 设备 性 能 指标 达到 
完全 失效 病 值 。 
失效 阐 值 
所 A eis B | C 
- : | 
wh! 
! 1 本 
图 1 典型 设备 退化 曲线 


Fig.1 Typical equipment degradation curve 
在 现 有 大 多 数 研究 中 , 国内 外 学 者 们 普遍 将 A、B、C 三 
阶段 退化 过 程 的 退化 因子 分 别 定义 为 常数 ， 线 性 函数 以 及 指 
数 函 数 。 对 于 指数 退化 阶段 ， 普 遍 将 相 邻 两 时 段 ! 与 1+At 的 
状态 内 转移 概率 表示 为 
w+D=aw(D4 ie(d23…N) 
其 中 乡 为 待 估计 参数 ， 为 定 值 。 基 于 Peng00 对 退化 因子 的 定 
义 ， 整 个 设备 加 速 退 化 为 一 递归 过 程 ， 对 于 离散 时 间 的 状态 
转移 ， 定 义 状态 退化 关系 式 为 1/ ， 递 归 闵 次 后 的 内 转移 概率 
值 为 1 ， 则 退化 函数 应 满足 : 
f°"(D) = ai(loa) (1) 
人 (2) 
式 (1) 表 示 刚 从 均匀 退化 状态 结束 时 设备 状态 内 转移 概 
率 ， 即 递归 的 初 值 ， 式 (2) 表 示 加 速 退化 过 程 中 加 速 的 概念 ， 
即 越 靠近 最 终 寿命 时 ,设备 停留 在 当前 状态 的 概率 下 降 越 快 。 
针对 加 速 退化 关系 式 的 定义 ， 在 分 析 指 数 递归 关系 式 结 
构 的 基础 上 ， 提 出 了 一 种 新 的 似 究 关系 式 加速 退 化 关系 : 
Qi(t+1)=G") 一 1 ie(l,2,3...N) (3) 
其 中 乡 为 待 估计 参数 ， 可 由 历史 数据 估计 得 出 ， 由 于 状态 转 
移 概率 os) 恒 小 于 1， 乡 为 1， 则 需要 减 去 边界 值得 到 适应 
后 的 状态 转移 概率 。 
对 于 选取 满足 递归 加 速 退 化 要 求 的 相 适 应 参数 进行 指数 
关系 与 似 震 关系 比较 ， 发 现 合理 参数 下 的 似 寡 关系 与 指数 关 
系 有 着 相似 的 递归 加 速 退 化 效果 ， 如 图 2 所 示 。 


姨 


0 - : 
12 03456 7 8 9 10 1 次 


(b) (d) 


图 2 关系 式 加 速 退化 递归 示意 图 


Fig.2 Schematic diagram of relational accelerated degradation recursion 
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其 中 子 图 (a)(b) 分 别 为 满足 加 速 退化 过 程 且 在 合理 参数 


刘 文 溢 ， 等 : 


训 


基于 改进 退化 隐 马 尔 科 夫 模型 的 设备 健康 诊断 与 寿命 预测 研 


下 的 指数 关系 及 似 需 关系 退化 示意 图 ， 粗 实 线 为 关系 式 的 连 
续 函 数 图 ， 粗 虚线 为 对 应 关系 式 的 反 函数 图 ， 点 线 为 递归 过 
程 中 状态 转移 概率 ,由 图 可 知 两 种 关系 式 具有 一 定 程度 的 作 
用 相似 性 。 子 图 (c)(d) 分 别 为 当前 示例 合理 参数 下 的 指数 关系 
退化 过 程 以 及 似 昧 关系 退化 过 程 ， 实 线 表示 各 个 单 时 点 状态 
概率 ， 虚 线 表 示 当 前 关系 式 下 概率 的 部 分 退化 过 程 。 由 图 看 
出 ， 似 究 关 系 在 表达 加 速 退化 过 程 中 ， 程 度 优 于 指数 关系 。 
2.2 退化 转移 概率 推算 

a) 平稳 退化 阶段 ， 相 邻 两 时 点 上 与 4+e 设备 保持 当前 状 
态 概 率 的 变化 为 定 值 ， 即 


oud+D-oai(D = 页 ie 人 2.3…N) 


2 


其 中 6<0， 保 持 当 前 状态 概率 @@ 减 小 ， 则 由 当前 状态 转移 
至 其 他 状态 的 概率 会 随 之 增 大 ， 且 根据 比例 分 配给 各 


(DD) CO=i+li+2…N)。 则 平稳 退化 状态 下 一 时 刻 状 态 转 移 概 率 为 
Qi(t+1)= 0;(1) + 
EAGAN 


oy (t+) =0()— (4) 
bp oD 
则 整个 平稳 退化 过 程 某 一 时 点 的 状态 转移 概率 可 由 初始 
时 刻 递 推 得 到 : 
Qi(kat) = (0)+ kG 
oD D+ (5) 
> oj(0) 


b) 均匀 退化 阶段 ， 相 邻 两 时 段 : 与 上 +Af 设备 保持 当前 
状态 的 概率 的 变化 为 ! 时 段 概率 的 线性 值 ， 即 
Git +D) -Q(t)=6, ilt) ie(2.3…N) 
其 中 <0，a@ 减 少 的 值 同样 分 配给 00 : 
ad+D=( 办 +Dai(D) 
Qj (Do 


> Oy (人 


jtl 


若 定 义 平稳 退化 阶段 结束 时 的 时 点 为 ， 则 处 于 均匀 退 
化 时 的 某 一 时 点 的 状态 转移 概率 可 描述 为 
Ota + kat) = (ph + Dt a(t) 


Gy (ta tka) = 0 (ta )+ Dh + — Do lie) 
Ya, (ta) 


c) 对 于 本 文 提 出 的 基于 窜 函 数 的 加 速 退 化 阶段 , 相 邻 两 
时 点 上 与 1+Af 设备 保持 当前 状态 的 概率 的 变化 为 1 的 似 窜 
关系 


Q(t+l)= 0(D)+ 


(6) 


(7) 


Qlt+D)=@" -1 iel(,2,3...N) 
其 中 >0,a) 减 少 的 值 分 配给 0 则 下 一 时 刻 的 状态 转移 


Calt+D)=p"" —1 
Oy (Db™™ -1- Ci())) 


Q(t+D)=0(D+ (8) 
2 %) 
同样 定义 均匀 退化 阶段 结束 时 点 为 刀 由 于 似 震 关系 特 


殊 的 递 推 形式 , 引用 递归 关系 符号 f*(W, 其 路 表示 递归 次 

数 ，X 表 示 递 归 概 率 初 值 ， 则 在 此 加 速 退 化 阶段 某 一 时 点 状 
Q(t + kat) = f+ (Q(t2) 

Q(t)(f (i(ta) -Gilt2)) 


Do 
过 设备 历史 数据 训练 后 得 到 初始 状态 转移 概率 矩 队 


Q(t + kat)= Gy (ta ) 十 


让 


[1 
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A ， 因 设备 的 运转 退化 过 程 在 无 人 工 干 预 的 自然 状态 下 是 不 
可 道 的 过 程 ， 即 设备 状态 只 由 停留 在 当前 状态 及 由 当前 状态 
退化 到 更 差 状 态 ， 因 此 初始 状态 转移 概率 矩阵 4 可 描述 为 
Qn ap i Qn Qn 
0 ao 1 Qn Co2N 
A=|: : | : 
0 0 ”CN_IN-L CN_IN 
0 0 .0 1 
结合 上 述 三 种 退化 递 推 式 及 转移 概率 矩阵 则 可 推算 出 
个 不 同 退 化 阶段 中 经 历 kat 后 的 状态 转移 概率 矩阵 为 
GOAG Or LDD 0)+ Sr OA 
> aij(0) > ou(0) 
4 -| 0 GO + + 
0 (10) 
0 0 : 
0 0 Re 1 
fa 
(0<k< Ee) 
a (b+) 四 -Do a (ph 4D -Da 
(办 +J) “on" ot+ E CO 十 
Zo; < 
”2 | 
0 0 
0 0 了 
( 亿 <K< 包 ) 
Ax Ax 
So) a + (G2) 0) 辣 a + ET (0) 8) 
2 之 
4 -| 。 a 加 CC 一) 
2 (12) 
0 0 : 
0 0 1 
(时 < 月 
上 述 式 (10) 表 示 从 初始 时 点 至 所 时 点 设备 系统 处 于 平稳 
退化 阶段 ， 经 过 kat (0<k< 章 ) 时 的 状态 转移 概率 矩阵 ， 描 述 
为 A 。 将 当 上 = 二 时 的 状态 转移 概率 矩阵 记 作 信 ， 描 述 为 


A,= (ww, 


阵 的 第 i 行 ? 
备 处 于 均匀 


其 中 史 为 平稳 退化 阶段 结束 时 状态 转移 概率 矩 
第 列 的 元 素 。 则 根据 描述 ， 可 得 式 (1D) 所 示 设 


退化 状态 , 上 且 经 历 ui gm 时 的 状态 转移 概 


率 和 矩阵 。 同 样 的 , 将 当 k= 过 时 的 状态 转移 概率 短 阵 记 作 A,， 


描述 为 4, =(@y)ww， 其 中 为 均匀 退化 阶段 结束 时 状态 转移 


概率 矩阵 第 


处 于 似 究 关系 加 速 退化 状态 并 经 历 


i 行 ， 第 j 列 的 元 素 ， [ 


3 ”基于 改进 遗传 算法 的 参数 估计 
对 于 极 大 似 然 估计 常用 的 EM 算法 容易 陷入 局 部 最 优 的 


缺陷 ， 


本 文采 用 了 全 局 搜索 能 力 更 强 的 遗传 算法 进 
DGHMM 的 参数 估计 ， 


T 


建 并 对 数 似 然 函数 
log(L(4, ,8,$)) =10g(Pr(O, 0..…0, [4,8,t, A)) 


此 得 式 (12) 所 示 系 统 设备 
万 了 kat 人 < 的 状态 转移 


交 


行 


202009.00081v1 
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将 每 组 用 以 迭代 的 参数 组 计算 而 出 的 似 然 值 作为 单个 解 
决 方案 的 适应 度 ， 则 每 一 次 种 群 迭代 的 最 优 个 体 便 可 表示 为 
4 人 =arg Imax log(L(4, $6.,) 


= arg max log(Pr(O, ,07 | ,6, t,t) 
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是 交换 随机 两 条 染色 体 各 自 一 部 分 位 点 ， 并 对 产生 的 新 染色 
体 进行 自 适 应 归 一 化 ， 交 叉 方 式 1、2 分 别 如 图 4、5 所 示 。 


el |el el je jej lelle 

el le el el je je je 

X lol le le le le lel le 

| 有 有 上 有 有 有 | [leHeHeHe He eel| 


在 经 过 若干 次 迭代 后 发 现 CO4., 欠 内, 办) 值 变化 不 大 或 达到 
预定 迭代 次 数 后 ， 便 输出 目前 较 优 的 参数 组 4 办 办 ,办 。 
3.1 改进 编码 

本 文采 用 混合 编码 策略 : 对 于 转移 概率 矩阵 ANxN) 以 及 
观测 抢 阵 BVNxM) ， 首 先 对 和 矩阵 进行 打 平 ， 由 于 转移 概率 拖 


| py 
\ VvV 


阵 自身 的 特性 ， 则 其 矩阵 中 实际 有 价值 的 数值 仅 有 


人 +N? 个。 考虑 到 概率 值 均 小 于 1， 且 大 部 分 为 多 位 数 的 
2 

小 数 ， 采 用 十 进 制 浮 点 数 编码 ， 且 保留 三 位 小 数 。 打 平 后 的 

转移 概率 年 阵 去 零 去 一 后 便 是 一 个 参数 组 中 转移 概率 染色 体 ， 

打 平 后 的 观测 概率 矩阵 便 是 参数 组 中 观测 染色 体 。 初 始 状态 

分 布 丈 同样 使 用 M 个 位 点 的 十 进 制 浮 点 数 编码 染色 体 表示 

(保留 三 位 小 数 )。 
对 于 入 ,为 ,考虑 到 整个 模型 对 退化 因子 极为 敏感 ,采用 


二 进 制 的 整数 编码 且 保 留 五 位 小 数 。 以 伪 18 个 位 点 的 二 进 
制 整 数 染色 体 对 54,8, 的 绝对 值 进行 转换 ， 其 中 第 一 个 位 点 
表示 参数 整数 部 分 ， 且 该 位 点 不 参与 二 进 制 与 十 进 制 的 相互 
转换 。 相 应 地 ,余下 的 十 七 位 在 解码 过 程 中 需要 缩小 十 万 倍 ， 
然后 加 上 首位 编码 数值 ， 便 得 到 解码 后 的 参数 。 各 变量 编码 
后 的 染色 体 如 图 3 所 示 。 

4 三 参 的 编码 中 ， 染 色 体 一 个 位 点 代表 了 和 矩阵 中 的 一 个 
概率 值 ， 每 个 位 点 由 四 个 部 分 组 成 。 灰 色 部 分 为 矩阵 中 单个 
值 的 整数 部 分 ， 白 色 部 分 分 别 为 该 值 的 十 分 位 、 百 分 位 以 及 
干 分 位 的 值 。 在 整个 改进 遗传 算法 流程 中 ， 单 个 参数 组 里 的 


参数 总 是 同时 参与 计算 , 一 个 参数 组 即 为 一 个 完整 的 解决 方案 。 
N2+N-2 
二 
i "EE 
NxM 
4 be bi J > 各 
Bp bm -- bo | 一 | 
百 分 | 
B= S we 
bn bi 加 | 至 执 基 的 做 村 世 
bn bm bm N 
A 
于 分 位 
而 分 位 
X= A 和 = 一 全 4 
“中 一 一 一 二 一 二 


整数 部 分 , 不 4 4 
/三 参与 进 制 转 换 。 17 位 二 进 抽 


4%% 一 人 十 TIITTTTTIT 


《保留 五 位 小 数 ) 


图 3 编码 示意 图 
Fig.3 Coding diagram 


3.2 ”改进 交叉 
对 于 上 述 改进 过 后 的 混合 编码 策略 ， 相 应 地 对 交叉 与 变 
异 方式 进行 改进 。 针 对 编码 方式 相似 的 多 三 参数 ， 设 计 两 种 


el lel|lellellellel le 
el lel|lellel le lel je | 
国 el le le jej le le je 
回国 本 = 天 古国 |eHjeHeHeHeHeHel| 
图 4 交叉 方式 1 
Fig.4 Crossover mode 1 
elj le le ljel lel lel le 
el je ljel ljel ljel ljel je 
图 DDNGNOGNOGUONC 
onggugg ee 加 国 国 时 国 国 国 
J 
el |el je je el je le 
el |el je je el lel je 
el |el je je el lel je 
HH HeHeHeHe He 王 e 有 二 
图 5 交叉 方式 2 


Fig.5 Crossover mode 2 
于 染色 体 是 概率 矩阵 打 平 后 的 产物 ， 珑 阵 同 行 元 素 和 
为 1 决定 了 自 适 应 归 一 化 需要 满足 一 定 范围 限制 ， 
适应 归 一 化 规则 如 表 1 所 示 ， 染 色 体 位 点 从 1 开始 计数 。 
表 1 自 适 应 归 一 规则 


Tab. 1 Adaptive normalization rule 


归 一 化 
染色 体 类 别 
A 组 别 归 一 范围 
1 1 一 N 
2 N+1—2N:1 
转移 染色 体 N*N) 
oe N+N-4 N+N-2 
2 2 
1 1 一 M 
四 M+1— 2M 
观测 染色 体 (N*M) 
N (N-DM+1—NM 
初始 染色 体 (N) 1 1 一 N 


4 三 参 的 变异 以 及 退化 因子 的 交叉 变异 均 同 常规 GA 算 
法 ， 变 异 操作 后 的 归 一 化 同样 遵从 表 1 规则 。 

改进 后 的 交叉 操作 ， 利 用 两 种 不 同 的 交叉 方式 ， 通 过 随 
机 选择 交叉 方式 ， 极 大 程度 地 增加 了 种 群 的 多 样 性 ， 增 强 了 
算法 的 全 局 搜索 能 力 ， 同 时 在 一 定 程度 上 保留 了 算 子 精度 ， 
有 利于 产生 全 局 最 优 参数 。 
在 算法 完整 运行 后 ， 使 得 似 然 函 数 取信 最 大 的 参数 组 
办 发 便 是 历史 数据 条 件 下 ，HMM 模型 最 优 的 估计 参数 。 
4 ”基于 DGHMM 的 设备 剩余 寿命 预测 
4.1 贪 林 近似 法 

对 于 常规 HMM 模型 的 预测 问题 ， 一 般 存在 近似 算法 以 
及 Viterbi 算法 。 现 阶段 研究 中 ， 普 遍 在 得 到 HMM 模型 参数 
以 及 观测 序列 后 运用 后 者 计算 此 观测 序列 下 概率 最 大 的 状态 


不 同 的 交叉 方式 1 与 交叉 方式 2。 交 义 方式 1 是 交换 随机 两 
个 染色 体 相同 位 点 上 千 、 百 分 位 数 和 十 分 、 整 数位 数 ， 并 对 
产生 的 两 条 新 的 染色 体 各 自 进行 自 适应 归 一 化 ; 交叉 方式 2 


序列 ， 即 
T=argmax Pr(@,0,…0 |4,7) 
实际 应 用 中 ， 在 需要 预测 设备 系统 剩余 寿命 时 ， 往 往 设 
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备 还 处 于 服役 中 ， 即 此 时 是 得 不 到 设备 从 服役 初期 开始 至 失 
效 的 整个 过 程 的 观测 值 。 为 了 解决 这 一 问题 ， 研 究 者 们 一 般 


be 


等 : 基于 改进 退化 隐 马 尔 科 夫 模型 的 设备 健康 诊断 与 寿命 预测 研究 


ww， 此 时 的 状态 转移 概率 久 
忆 子 及 初始 状态 转移 概率 矩阵 
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E 阵 中 第 1 行 的 数据 可 由 上 述 退 化 


E 算 而 知 ， 则 设备 的 剩余 寿命 


用 同 种 设备 历史 数据 加 以 代 蔡 并 建立 模型 ， 输 入 服役 设备 的 。” ”可 表示 为 
某 一 时 点 之 前 的 观测 值 序列 ， 预 测 其 隐 状 态 ， 再 通过 役 龄 时 Ri 5 
间 分 布 来 确定 剩余 寿命 ， 但 这 与 HMM 模型 必须 服从 指数 分 ™ 六 or 名 (14) 
布 的 局 限 性 冲突 ， 且 设备 之 间 存 在 的 固有 差异 以 及 设备 状态 全 
的 随机 性 ， 使 得 预测 在 一 定 程度 上 损失 设备 适应 性 。 式 (14) 中 乙 为 历史 数据 中 设备 处 于 7 状态 的 寿命 ,， as” 为 
寻 此 本 文 在 近似 算法 以 及 Viterbi 算法 基础 上 提出 一 种 ”当前 状态 i 最 后 一 个 节点 时 的 由 i 状态 转向 i 状态 的 概率 值 ， 
预测 剩余 状态 的 新 方法 (下 称 贪 禁 近似 法 )。 ay” 为 当前 状态 i 最 后 一 个 节点 时 的 由 i 状态 转向 /状态 的 概 
贪 禁 近 似 法 将 近似 算法 寻求 状态 的 思想 ， 运 用 到 寻找 剩 。 率 值 。 即 设备 总 剩余 寿命 等 于 当前 状态 寿命 余 量 以 及 后 续 所 
余 观测 上 ， 即 在 已 知 HMM 参数 4 的 前 提 下 ， 寻 求 下 一 时 刻 有 状态 的 寿命 之 和 。 
最 有 可 能 出 现 的 观测 值 。 ‘er 
定义 设备 在 第 1 时刻 产 生 观 测 & 的 前 提 下 ， 第 !+1 时 刻 5 算 例 分 析 
产生 观测 9 的 概率 为 .0=Pl0,=90.10=0,D (k,g =12…M)， 通过 美国 卡特 彼 勒 公司 液压 泵 的 设备 健康 诊断 与 寿命 预 
则 Y 的 数学 描述 为 测 实 例 来 验证 评价 本 文 提 出 的 模型 与 方法 。 实 验 室 中 设备 的 
振动 信号 是 由 安装 在 与 液压 泵 旋转 轴 平 行 位 置 的 液压 加 速 计 


(k,g =1,2:….M) 


OI A 


(13) 
he 


其 中 为 :时 刻 设备 系统 的 状态 转移 概率 ， 其 如 上 述 所 说 符 
三 阶段 退化 代入 不 同 的 k,g 值 便 能 计算 出 从 1 时 刻 的 大 观 


测 转移 到 ! 十 1 时 刻 的 8& 观测 的 概率 , 由 此 得 出 间接 观测 转移 
Vi WV, Vu Vi 
WW Vs Vu Van 
V'= 
Vi Vi Vai Vi 
Va Vu i Viam Vum 


于 状态 转移 矩阵 4 随 着 监测 时 点 逐渐 退化 , 则 WV 也 同 
样 是 随 着 时 间 变 化 的 。 对 于 某 一 产生 观测 8(=42…M) 的 时 
点 !， 其 下 一 时 点 1 十 1 最 有 可 能 出 现 的 观测 值 概率 为 


P(t+D)= maxv(i=1,2.…M) 
<j<, 


下 一 时 点 最 优 可 能 观测 则 为 


0 1()) =arg max vi (i=1,2.…M) 


I+ 


ps 列 为 前 t 个 实际 观测 值 Cu =(0,,0,…0,) 


“0 0) 。 


I 步骤 如 下 ; 

a) 将 Ow 代入 Viterbi 算法 求 得 最 大 概率 状态 路 径 Roadw ， 
并 记录 最 后 一 个 状态 值 为 * 。 

b) 运用 近似 法 求 得 下 一 时 点 最 有 可 能 出 现 的 观测 2 
将 其 拼接 到 0 之 后 形成 Cs 。 

c) 将 Oi 代入 Viterbi 算法 求 得 新 的 最 大 概率 状态 路 径 
Road,。, ， 并 记录 最 后 一 个 状态 值 为 * 。 

qd) 判断 * 与， 如 果 s=s*, 则 重复 步 马 b)~d); 如 果 ss 
则 至 步 又 e)。 

e) 最 初 的 0w 与 最 终 的 2 相 减 便 是 预测 的 剩余 观测 ， 
输出 剩余 观测 。 

在 求 得 剩余 观测 之 后 ， 运 用 标准 Viterbi 算法 , 代入 加 入 
剩余 观测 后 的 观测 序列 及 模型 参数 ， 寻 求 当 前 子 状态 下 概率 
最 大 剩余 状态 路 径 fw 。 

4.2 设备 剩余 寿命 预测 

对 于 4.1 节 求 得 的 最 大 剩余 状态 路 径 fn 假设 当前 设备 

处 于 状态 i， 则 在 Zs 中 找到 处 于 状态 i 的 最 后 一 个 时 间 节 点 


收集 。 在 应 用 实例 中 , 分 别 对 液压 泵 充 入 20、40、60 与 80mg 
10min 采集 一 个 时 间 大 约 为 1min 的 振动 信 


的 微 尘 ， 


实验 过 程 中 ， 液 压 
Cont2 以 及 Cont3， 


Cont1 、 


并 每 隔 
号 样本 。 随后 使 用 


10db 的 小 波 将 振动 信 
降 维 后 的 小 波 系数 作为 DGHMM 的 输入 特征 序列 向 量 。 


号 分 为 五 层 ， 


将 经 过 


整个 


泵 的 状态 可 分 为 


种， 分 别 为 Baseline、 


其 中 Cont3 状态 为 设备 的 彻底 失 


整个 实验 分 析 平 台 为 Python3， 平 台 运行 环境 为 Windows 10。 
5.1 DGHMM 数据 准备 


利用 来 自 36 个 传感器 (Sensor) 的 振动 信号 和 对 


液压 泵 进行 健康 状态 诊断 以 及 剩余 寿命 预测 。 


部 分 经 过 小 波 


变换 后 的 振动 数据 (Pump6) 如 表 2 所 示 。 
表 2 Pump6 部 分 小 波 变 换 数 据 


Tab.2 Partial wavelet transform data of pump6 


时 Sen Sen Sen Sen Sen 
点 1 7 15 23 32 

1 2.62 19.53 5.66 0.06 0.96 
3 2.65 20.24 6.10 0.06 1.04 
5 2.39 17.45 4.92 0.06 0.87 
起 2.07 16.15 3.98 0.06 0.75 
9 2.25 18.53 4.51 0.06 0.81 
11 2.24 21.44 4.38 0.06 0.87 
13 2.41 5.37 6.84 0.08 1.66 
15 2.58 6.16 7.96 0.09 1.84 
17 2.53 6.07 7.63 0.09 1.84 
19 2.44 6.24 7.71 0.10 1.84 
21 20.66 9.68 8.42 0.11 2.17 
23 5.06 8.75 7.20 0.10 1.89 
25 7.53 8.82 8.15 0.10 2.00 
27 5.74 8.37 7.19 0.10 1.89 
29 3.46 8.25 7.16 0.10 1.92 
31 41.78 7.98 7.12 0.15 1.86 
33 31.08 7.90 6.98 0.14 1.84 
35 60.15 7.99 7.03 0.15 1.85 
37 了 7.27 7.80 7.12 0.17 1.84 

5.2 实验 参数 估计 
本 文 运用 改进 遗传 算法 同时 对 DGHMM 与 时 变 HMM 进 


行 参数 估计 。 表 3 与 表 4 分 别 为 泵 
效 临界 状态 转移 概率 矩 
为 0.7， 变 异 概率 为 0. 


阵 。 算 法 中 种 群 数量 为 30， 
3。 考 虑 到 不 同 参数 组 下 实验 数据 出 现 


6 估计 参数 组 中 初始 与 失 
交叉 概率 


的 概率 似 然 值 可 能 


存在 数量 级 差异 较 大 的 情况 ， 通 过 减少 种 


群 数 量 、 增 加 进 代 次 数 来 增强 选择 力度 ， 相 应 地 通过 增加 变 
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异 概率 来 保证 每 代 种 群 多 样 性 。 表 6 泵 6 与 泵 82 在 各 模型 下 路 径 概率 值 
DGHMM 与 时 变 HMM 在 改进 GA 算法 整个 流程 中 的 概 Tab.6 The path probabilities ofpump 6 and pump 82 under each model 
率 似 然 值 迭 代 曲 线 如 图 6 所 示 (Pump6)。 从 参数 估计 的 结 类 别 Pump6 Pump82 
来 看 , 实验 数据 在 DGHMM 的 较 优 参数 组 下 出 现 的 似 然 值 大 HMM 7.32E-11 8.16E-11 
于 在 时 变 HMM 的 较 优 参数 组 下 出 现 的 似 然 值 。 这 一 结果 在 上 时 变 HMM 5.39E-13 1.03E-13 
一 定 程度 上 可 反映 DGHMM 较 常 规 基于 指数 式 加 速 退化 的 DGHMM 1.36E-13 1.77E-12 
时 变 HMM 而 言 ， 能 更 好 地 描绘 实验 数据 。 5.4 剩余 寿命 预测 
表 3 Pump6 初始 状态 转移 矩阵 运用 提出 的 贪 禁 近 似 算 法 以 及 常规 Viterbi 算法 的 设备 
Tab.3 Initial state transition matrix of pump6 剩余 寿命 预测 方法 ， 动 态 地 求 得 Pump6 与 Pump82 在 
state baseline contl cont2 cont3 DGHMM 下 的 剩余 寿命 ,结果 由 图 7 所 示 。 从 图 中 可 以 看 出 
baseline 0.924 0.064 0.007 0.005 预测 而 出 的 剩余 寿命 值 与 实际 剩余 寿命 的 误差 处 于 相对 合理 
contl 0.000 0.783 0.128 0.089 范围 内 ， 且 当 设 备 处 于 失效 状态 cont3 时 ， 设 备 停 止 运作 ， 
cont2 0.000 0.000 0.807 0.193 退化 相应 停止 ， 模 型 也 相应 地 成 功 识 别 并 判定 设备 剩余 寿命 
cont3 0.000 0.000 0.000 1.000 为 0， 虽然 预测 值 一 定 程度 上 在 Baseline 与 contl 状态 转换 
表 4 Pump6 失效 临界 状态 转移 矩阵 处 存在 “滞留 现象 "， 但 整体 而 言 ， 基 于 DGHMM 的 贪 禁 近 似 
Tab.4 Failure critical state transition matrix ofpump6 法 对 于 Pump6 与 Pump82 剩余 寿命 的 预测 是 有 效 可 行 的 。 
state baseline contl cont2 cont3 本 
baseline 0.426 0.483 0.053 0.038 | 
contl 0.000 0.361 0.377 0.262 200- 
TT cont2 0.000 0.000 0.403 0.597 县 Ei 
> cont3 0.000 0.000 0.000 1.000 证 


一 一 实际 了 LL， 一 一 预测 


到 


图 7 Pump6 与 pump82 剩余 寿命 预测 
Fig.7 RUL ofpump6 and pump82 


6 ”结束 语 


本 文 提 出 了 一 种 以 新 加 速 退化 为 核心 的 设备 健 


对 数 似 然 度 


康 诊断 与 


剩余 寿命 预测 方法 。 提 出 的 似 寡 关 系 更 保守 ， 较 同 为 退化 医 

1 而 am a 2 子 的 指数 式 关系 而 言 更 稳定 ， 且 能 更 好 描述 设备 性 能 随 役 龄 

进化 代数 曾 加 的 逐步 加 速 退化 的 过 程 ,对 于 设备 的 状态 识别 也 更 精确 ; 

一 一 DOHWM ~ 时 变 HWN 车 于 改进 遗传 算法 的 参数 估计 方法 ， 也 取得 了 相对 较 好 的 效 

图 6 Pump6 似 然 闪 代 对 比 图 果 ， 提 出 的 贪 禁 近似 法 一 定 程度 上 打破 了 常规 隐 马 尔 可 夫 模 

Fig.6 Likelihood iteration comparison diagram of pump6 型 时 间 上 必须 服从 指数 分 布 的 束缚 ， 能 够 近似 地 对 设备 的 剩 

5.3 ”健康 状态 诊断 余 寿 命 进 行 预测 。 因 此 基于 改进 退化 隐 马 尔 可 夫 模型 的 设备 
分 别 将 泵 6 与 泵 82 的 实验 数据 向 量化 后 输入 HMM、 时 ”健康 诊断 与 剩余 寿命 预测 方法 在 描绘 设备 退化 、 设 备 状 态 诊 


ea 


变 HMM 以 及 DGHMM 模型 ， 运用 Viterbi 算法 分 别 求 得 断 准 确 率 方面 更 加 有 效 ， 在 剩余 寿命 预测 上 亦 为 可 行 。 
在 各 自 实验 观测 数据 和 各 自 较 优 参数 组 下 设备 的 最 大 概率 本 文 仅 将 似 容 关系 应 用 到 了 HMM 模型 并 验证 了 结果 ， 
态 路 径 并 加 以 验证 评价 ， 结 果 如 表 5 所 示 ， 泵 6 与 泵 82 实 未 来 应 将 其 应 用 到 HMM 的 拓展 一 -HSMM 中 , 进一步 对 似 


> 


验 数 据 在 不 同 模型 下 的 最 大 概率 状态 路 径 概率 值 如 表 6 所 示 。 窜 关 系 进行 验证 。 为 保证 寻 优 效果 ， 文 中 的 改进 遗传 算法 连 
表 5 泵 6 与 泵 82 健康 状态 识别 结果 代 次 数 较 多 ， 时 间 复 杂 度 较 高 ， 未 来 也 应 尝试 改善 ， 且 在 贪 
Tab.5 Health states identification results of pump6 and pump82 焚 近 似 法 进行 过 程 中 若 出 现 无 限 循环 情况 ， 还 应 设置 相应 的 
Dili pump82 和 迭代 次 数 限制 才能 完成 算法 相应 步骤。 
状态 base contl cont2 cont3 准确 率 状态 base contl cont2 cont3 准确 率 
base 12 0 0 0 100% base 9 0 0 0 100% 参考 文献 : 
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oe cont2 0 3 8 0 73% con2 0 4 6 0 60% 计算 机 研究 与 发 展 , 2004, (08): 1421-1429. (Liu Qun, Zhang Huaping, 
cont3 0 0 0 8 00% cont3 0 0 1 11 92% Yu Hongkui, et al. Chinese Lexical Analysis Using Cascaded Hidden 
base 12 0 0 0 00% base 9 0 0 0 100% Markov Model [J]. Journal of Computer Research and Development, 
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